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1. Introduction 



Non-hermitean random matrices have interesting and important apphcations in various 
physical problems. Examples for such applications are the phase diagram of QCD[4[Ej, 
nuclear decay and resonances in multichannel scattering in chaotic and disordered 
systems jHl Ej, neural networks [Zj, and asymmetric random hopping matrices - the 
Hatano-Nelson model |H] - with applications to pinning of magnetic flux lines in high 
temperature superconductors and also as a new tool for studying Anderson localization 
in disordered conductors. Studying Anderson localization by means of non-hermitean 



Non-Hermitean RMT: Planar Diagrams and the "Single-Ring" Theorem 



2 



matrices is also one of the reasons for the very recent growing interest in studying the 
statistics of resonance widths and delay times in random media OCD] . Our final example 
is diffusion in random velocity fields which gives rise to random Fokker-Planck 
operators. For a nice recent overview of random matrix theory and its applications, 
including non-hermitean matrices, see jT2]. 

1.1. The Difficulty Posed by Non-Hermitean Matrix Ensembles 

Let us recall first a few basic facts about hermitean random matrices. A basic tool 
in studying hermitean random matrices (henceforth all matrices will be taken to be 
N X N with N tending to infinity unless otherwise specified) is the Green's function 
defined by 



where (■ ■ ■) denotes averaging over the random distribution from which the matrices 
are drawn. Diagonalizing by a unitary transformation we have 



where the real numbers {A^} are the eigenvalues of 0. Thus, G{z) is a meromorphic 
function with poles at the eigenvalues of 0. In the large N limit, the poles merge into a 
cut (or several cuts) on the real axis|. Powerful theorems from the theory of complex- 
analytic functions can then be brought to bear to the problem of determining G{z) 
All of this is well-known. 

In contrast, for a non-hermitean matrix, the eigenvalues invade the complex 
plane. For example, for non-hermitean matrices generated with the probability 
— ;|e~^*'''^^'^, Ginibre determined long ago ^3] that in the large-A^ limit, the 
eigenvalues are uniformly distributed over a disk of radius unity in the complex plane. 
The Green's function corresponding to this ensemble is 



and it is obviously not analytic inside the unit disk. In general, the eigenvalues of a 
non-hermitean matrix taken from some probability ensemble, will occupy, on average, 
a two-dimensional domain T> of the complex plane, where the Green's function is non- 
analytic. Thus, we lose the powerful aid of analytic function theory in analyzing these 
ensembles. 

In the theory of hermitean random matrices, a powerful method consists of the 
Feynman diagrammatic expansion ,15. in which one expands G{z) = 

I More precisely, this is really a statement on the eigenvalues of an individual realization of cj) in the 
large N limit (assuming that all eigenvalues are typically confined in a finite segment, independent of 
N). The averaging process smudges the eigenvalues into cuts already at finite N. 




(1) 





z\ < 1 



z\ > 1 



(2) 
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a series in 1/z, the " bare quark propagator." This method is clearly no longer available 
in studying non-hermitean random matrices. The knowledge of G{z) as a series in 
1/z can no longer tell us anything about the behavior G{z) for small z, as the simple 
example Q above makes clear. The eigenvalues fill a two-dimensional region rather 
than a one-dimensional region, as is the case for hermitean matrices. 

A renormalization group inspired method used in JH] implicitly involves an 
expansion in 1/z and is thus also not available for dealing with non-hermitean matrices 
without suitable further developments. 



1.2. Some Basic Formalism and Notations 

I shall now derive a couple of standard formulas which play a central role in studying 
the eigenvalue distribution of random non-hermitean matrices. The derivation is 
straightforward (see, e.g., [7j and Section 2 in P, and references therein), and is based 
on the identities§ d {1/z*) = d* {1/z) = Tx5{x)5{y) and dd* \og{zz*) = '7i6{x)6{y) . 

We are now ready to deal with the average density of eigenvalues of a non-hermitean 
matrix 0. We diagonalize the matrix by a similarity transformation || 

= S-^AS (3) 

where A denotes a diagonal matrix with elements Aj, i = 1, ...,N. Taking the hermitean 
conjugate we have (j)^ = S'^A*S^^'^. By definition, the density of eigenvalues is 

p{x, 2/) = E - ^^«) - ^^^)) • (4) 

i 

From the identities mentioned above, we obtain that d* tr {z — A)~^ = d* tr {z — 
= Tr^.5{x - ^Xi) 5{y - ^\,). Thus, using the identity det{z - A){z* - A*) = 
det{z — (j)) {z* — 0^) we may express p in terms of (p, and in a manner which is manifestly 
symmetric in z and z*, as 

p{x,y) = ^dd* (^tr log {z - <jy){z* - 0^)) . (5) 

This expression involves the logarithm inside the average, as well as two derivatives. 
Thus, unless one has a simple way of calculating the average in ^ {e.g., the replica 
method in the case of Gaussian ensembles [7j), (0) is not a practical expression for 
p{x,y). In many cases it is easier to calculate the Green's function 

(with \i = Xi + iyi) than to do the average (;^tr log {z — (f>){z* — 0^)) . Then, from (jH) 
and from one of the identities mentioned above, we have 

p{x,y) = -d* G{z,z*), (7) 

71 

§ Our notations are standard: ioi z — x + iy we define d = = ^ (^-^ — i-§^ so that dz — 1. Clearly, 

dz* = 0. Similarly, we define d* = ^ = \{j^+ z^) , so that d* z* = 1. 

II Non-diagonalizable matrices constitute a set of measure zero in the ensemble. 
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which is a simpler representation of p{x,y) than (0). 

These two representations for p{x, y) and the relation between them can be 
interpreted by recognizing (jSl) as a two dimensional Poisson equation for the electrostatic 
potential — tr log {z — (j)){z* — 0"'")) created by the charge density p{x,y). This 
connection between eigenvalue distributions of complex matrices and two-dimensional 
electrostatics has been long known in the literature [7j. Continuing along this line, 
consider the Green's function G{z,z*) in ©• If we define the electric field E = 
(SftG, —QG) then from the definition of G{z, z*) we have 

E(x) = / rf^2/p(y)^f^ (8) 

and thus ((Tj) is simply the statement of Gauss' law for this electrostatic problem. 



1.2.1. A Quaternionic Interlude and Speculation I end this section with a speculative 
remark, made for the first time in pQ. As was pointed out earlier, in the theory of random 
hermitean matrices it is well known that it is much easier to work with G{z) = 
rather than with the density of eigenvalues p{fi) directly, since the power of analytic 
function can be brought to bear on G{z). Of course, one can go from p(x) to G{z) with 
the identity 

- = ^ Six) . (9) 

Here, for random non-hermitean matrices our formula (0) gives the density of 
eigenvalues p{x, y) directly, but this formula is awkward to work with because of the 
logarithm. Its direct analog in the theory of random hermitean matrices would be 
p{x) = -^{jjtr9{x — (f))), which would also be extremely awkward to work with. It 
would thus be desirable to define a quantity analogous to G{z) and develop a method 
for calculating it. Reasoning along these lines we are led to an attempt to write a 
function of a quaternionic variable, in the same way that going from p{x) to G{z) we 
went from a function of a real variable to a function of a complex variable. Indeed, the 
obvious analog of Q is easy to find, namely. 



1 

27 



1 



[z-je)\z-je\ ^ \z - je)\z - je\ 

27rS{x) S{y) . (10) 



Here z = x + iy is a complex number, e is a small positive real number and {1, i,j, k} 
is the standard basis of the quaternion algebra^. (The absolute value of a quaternion 
is defined by \a + ib + jc + kd\ = ^/a^~+l^~+~cF^^~cP .) 

% One may of course write an equivalent formula with j replaced everywhere by the third quaternion 
basis element k. Also, e may be taken complex. 
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Unfortunately, (fTUI) as it stands, is less useful than because it leads to the 
quaternionic Green's function G{q) {q being a quaternionic variable) 



which involves the absolute value operation explicitly, and thus cannot be written as a 
simple trace like (Q). Note, however, that if we split q quite generally into q = z + jw 
{z, w being ordinary complex variables), then G{z + jw) is manifestly non-analytic in z 
as w — > 0, even before taking the average. 

This was posed in [IJ as an interesting problem, namely to find a useful quaternionic 
generalization of (0) for random non-hermitean matrices, and to develop the analogue 
of the work of Brezin et al. ^Hj associated with it. Unfortunately, this problem is still 
open, at the moment of transcribing this talk. 

The rest of this paper is organized as follows: In Section 2 I present the Method 
of Hermitization [H El El ^] by which the problem of determining the eigenvalue 
density of random non-hermitean matrices can be reduced to the more familiar problem 
of determining the eigenvalue density of random hermitean matrices. In Section 3 
I explain how this method was applied in j2j to derive a master formula for the 
density of eigenvalues of non-hermitean random matrices taken from a large class 
of rotationally invariant non-gaussian probability ensembles. This section contains a 
detailed pedagogical exposition of the diagrammatic method and its large- behavior, 
which is the basis for the discussion in that section. In Section 4 I present numerical 
simulations of the quartic ensemble, carried in jS], and demonstrate that they fit well 
with the theoretical predictions of the formalism developed in Section 3. In Section 5 I 
formulate and prove the "Single Ring" Theorem of [2], according to which, in the large- 
limit, the shape of the eigenvalue distribution associated with any of the ensembles 
studied in Section 3 is either a disk or an annulus. Finally, in Sections 6 and 7, following 
I present universal properties of the disk and annular phases and of the transition 
between them, as well as nice explicit expressions for the location of the boundaries of 
the eigenvalue distribution and its boundary values. Unfortunately, for lack of space, 
I cannot cover the interesting topic of adding random non-hermitean matrices. The 
interested reader could read about it in 

2. The Method of Hermitization: Reduction to Random Hermitean 



I have emphasized that a straightforward diagrammatic method is not allowed for non- 
hermitean matrices. Somewhat remarkably, it is possible to arrive at a diagrammatic 
method indirectly. I will show that the problem of determining the eigenvalue density 
of random non-hermitean matrices can be reduced to the problem of determining the 
eigenvalue density of random hermitean matrices, for which the diagrammatic method 
may be applied. This is achieved by the Method of Hermitization, as it was dubbed in 




(11) 



Matrices 
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Variants of this idea were presented independently in 013 HI]- (Here we follow the 
conventions of ^J.) 

We start with the representation p(a;, ?/) = ^dd* {jfti log {z — (j)){z*—(f)^)) (namely, 
Eq. (0).) By standard manipulations we observe that (tr^^^ log(2; — (p){z* — 0^)) = 
(tr^j^j logH) — inN, where H is the hermitean 2N x 2N matrix Hamiltonian 



H 



( 



\ 



\ 



Thus, if we can determine 

i^(^;^,^*) = ^(tr(,,) \og{n-E)) 

we can determine p{x,y). 

The matrix H is what is called a chiral matrix 

1 



(12) 



(13) 



anticommutes with the 'Tc;" matrix 



Due to its block structure, it 
, with the immediate consequence 



-1 

that its spectrum consists of pairs of eigenvalues with equal magnitudes and opposite 
signs. Namely, if ^ is an eigenvalue of H, so is — ^. 

Consider now the propagator matrix associated with H, namely, 

1 



g^^{T];z,z*) = ( 



) 



(14) 



where ?7 is a complex variable and the indices /i and u run over all possible 2N values. 
Here we followed a common practice and borrowed some terminology from quantum 
chromodynamics: we may consider (f),(j)^ as "gluons" (in zero space-time dimensions), 
which interact with a 2N dimensional multiplet of "quarks" ip^, with a complex mass 
matrix (the bare "inverse propagator") 



( 



\ 



V 



\ 



V J 



(15) 



(expressed in terms of its iV x blocks). Q^i, is thus the propagator of these quarks in 
the fluctuating gluon field. The crucial point is that since H is hermitean, Q^^^, can be 
determined by the usual methods of hermitean random matrix theory. In particular, as 
we already mentioned, the diagrammatic evaluation of ()14|) is essentially the expansion 
of Q^i, in powers of 1/r], with interaction vertices H. This is a well defined procedure for 
large rj, and it converges to a unique function which is analytic in the complex rj plane, 
except for the cut (or cuts) along the real axis which contain the eigenvalues of H. After 
summing this series (and thus determining Q^uiv'i -^j -2*) in closed form), we are allowed 
to set ?7 — > in ()14|) . Speaking colloquially, we may say that the crucial maneuver here 

^ The problem of determining the Green's function of chiral matrices such as H has been discussed by 
numerous authors ^I^J^]. (See also the paper by Feinberg and Zee cited in 116,1.) 
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is that while we cannot expand in powers of 2, we can arrange for z to "hitch a ride" 
with rj and at the end of the ride, throw rj away. 

We can now calculate p{x, y) by two different methods, each with its advantages. 
The first is to observe that 

/O ^\ 



1 

H 



J 



Thus, the quantity is simply the lower left block of Q^uiv = 0; z, z*). In other words, 
once we have Q^uiVi ^5 ^*) we can set rj to zero and use (jHl) and ((Zj) to write 



An alternative is to take the trace of ^^^^ 

1 , 1 



:(tr 



) 



G{0;z,z* 



(tr 



1 



(iV) ^2_(^*_0t)(^ 



(16) 



from which one can determine (jl3|) (in terms of a dispersion integral), and therefore 
determine p{x,y), as explained in [1^ in detail. I shall not pursue this possibility here. 
However, a few important comments relating to (jl7|) are in order. Observe from (jl7|) 
that 

^(0 - zO; z, z*) = ^(tr^^, 5 [^{z* - 0t)(^ _ 0))) (18) 

and thus counts the (average) number of zero-eigenvalues of the positive semi-definite 
hermitean matrix a/ {z* — (t)^){z — 0). Thus, if (on the average) has no eigenvalues 
equal to z, or in other words, if the density of eigenvalues of (p vanishes at z, then 

g{0-iO;z,z*) = 0, (19) 

independently of the large N limit. In particular, the boundaries of the eigenvalue 
distribution of in the large limit, are the curve (or curves) in the complex plane 
which separate regions where Q{0;z,z*) = from regions where Q{0;z,z*) 7^ 0. As a 
matter of fact, one can infer the location of these boundaries even without an explicit 
knowledge of ^(0; z, z*) and G{z, z*), by investigating the "gap equation" for ^(0; z, z*), 
namely, the trace of both sides of (jHl), and then setting r/ = p. The "gap equation" is 
a polynomial in Q with real coefficients (which depend on z explicitly as well as through 
G{z,z*), and the parameters which appear in the probability distribution for 0.) Due 
to the chiral nature of H, this polynomial contains a trivial factor of Q{0]z,z*) which 
we immediately factor out. Setting ^ = in the remaining factor we obtain an equation 
for the boundary. This is discussed in detail in p. 

It is worth emphasizing that the formalism developed here has nothing to do with 
large A^ as such. It is also totally independent of the form of the probability distribution 
P(0, 0^). The formalism is of course indifferent to the method one may choose to use 
to determine Q^y{rj; z, z*). 
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To summarize, we have obtained a formalism - the "Method of Hermitization" - 
for reducing the problem of dealing with random non-hermitean matrices to the well- 
studied problem of dealing with random hermitean matrices. Given the non-hermitean 
matrix 0, we study the hermitean matrix H instead. By whatever method one prefers, 
once one has determined the quark propagator Qf^i, (or its trace, the Green's function 
Q{ri] z, z*)), one can in principle obtain p{x,y) using (fTBj) (or using the other method 
mentioned following (fT7j) ). Whether that can be done in practice is of course another 
story. In the next section I shall present a large family of probability distributions 
P{4>, 0^) for which it is possible to compute Q^uijl'i ^i ^*) closed form (in the large- iV 
limit) by a diagrammatic method. Therefore, for those models, the algorithm discussed 
above can in fact be followed all the way to the end. 

3. Rotationally Invariant Non-Gaussian Ensembles 

The literature on random non-hermitean prior to j2] has focused exclusively on Gaussian 
randomness, a prime example being Ginibre's work J3] on the probability distribution 
P(0) = (1/Z)exp {—Nii(f)^(f)). In [2] it was shown that using the method of hermitization 
one can determine the density of eigenvalues of probability distribution of the form 



where V is an arbitrary polynomial of its argument. Indeed, by a simple trick, it was 
shown in [2] that one can obtain the desired density of eigenvalues with a minimal 
amount of work, by judiciously exploiting the existing literature on random hermitean 
matrices. In this section I explain how this can be done. 

In a certain sense, the work presented in pj may be thought of as the analog of 
the work of Brezin et al. for random hermitean matrices jT^j; they showed how the 
density of eigenvalues of hermitean matrices taken from the probability distribution 
P(0) = (1/Z)exp[— A^tr\^(0)] with V an arbitrary polynomial can be determined, and 
not just for the Gaussian case studied by Wigner and others, in which V = (l/2)tr0^. 
An important simplifying feature of the analysis in ^3] is that P{(f)) depends only on 
the eigenvalues of 0, and not on the unitary matrix that diagonalizes it. In contrast, 
the probability distribution (j20|) for non-hermitean matrices depends explicitly on the 
GL{N) matrix 5* used to diagonahze = S^^AS, and 5* does not decouple. Remarkably 
however, for the Gaussian P{<f)), Ginibre ^3] managed to integrate over S explicitly and 
derived an explicit expression for the probability distribution of the eigenvalues of 0. 
Unfortunately, it is not clear how to integrate over S and derive the expression for the 
eigenvalue probability distribution for non-Gaussian distributions of the form ()20|1 . In 
[2j this serious difficulty was circumvented by using the method of hermitization. 

Due to the symmetry of P{4>) under the transformation — > e*"0, the density 
of eigenvalues is obviously rotational invariant. It is useful to introduce to circularly 
invariant density as 



P(0) 



-NtTV{4>''(f>) 



(20) 



(21) 
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Rotational invariance thus leads to a simpler form of the defining formula (jH)) for G{z, z*), 
reproduced here as 



2^/ Pix',y') 



z — z' 



which reads 



7(r) = zG{z, z*) = / r'dr' p(r') , (22) 







whence 



p(r) = . (23) 
r ar 

Clearly, the quantity 7(r) is a positive monotonically increasing function, which satisfies 
the obvious "sum-rules" 

7(0) = and 7(00) = 1 . (24) 

In particular, observe that the first condition in ()24|) insures that no S{x)5{y) spike arises 
in p{x,y) when calculating it from ((7j) with G{z,z*) given by (j^ . as it should be. We 
are now ready to perform the diagrammatic expansion for probability distributions of 
the form 

3.1. The Diagrammatic Expansion and its Resummation 
Consider the matrix Green's function Q^i, in (fT^ and write it as 

where Qq^ is the inverse propagator defined in ()15|1 and V is the random chiral matrix 
(the "gluon field") 

/ 

(26) 

y 0t y 



V 



For large complex rj we may perfom the perturbative expansion of Q^j^y in powers of Qq 

Q,. = Qo,u + {{Qovg,)^;) + {{Q,vgovg,)^j) + ■■■ (27) 

which converges to a unique matrix function, analytic in the //-plane, cut along the 
appropriate segments on the real axis. This was our motivation to "hermitize" the 
problem to begin with. 

The diagrammatic method ^Hl provides a convenient and efficient way of 
representing the various terms which appear in the perturbation series ()27|). 

* Obviously, only terms with an even number of y's survive in H27f) for probability distributions of the 
form (EOl)- 
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Before sketching the Feynman rules for this expansion, some important symmetry 
considerations are in order. 

As was mentioned following Gfj,i, is essentially the propagator for the quarks 

ipf^ {fi = 1,2, ... , 2N) in the fluctuating gluon field V. Let us spht the 2N dimensional 
quark multiplet tp^ into two N dimensional "flavors", ip = {u,d). We see then that the 
"quark-quark-gluon" interaction is 



This interaction term, and the probability distribution ()2()|1 are manifestly invariant 
under the gauge transformation — > P4>Q\ (p^ Q(f)^P^, u Pu and d Qd, where 
P and Q are two independent unitary matrices, i.e., under the U{N) © U{N) subgroup 
of U{2N). These transformations may be written more compactly as 



We shall index quantities which transform under the first U{N) in (j^^ by Latin indices 
i,j,..., and quantities which transform under the second U{N) by indices from the 
beginning of the Greek alphabet a, (3, . . ..(Indices from the middle of the Greek alphabet, 
. . . are reserved for 2N dimensional objects.) Thus, for example, we have Ui, do, and 
the mixed object (pia . 

3.1.1. Feynman Rules and the Topological Nature of the Expansion Let us turn now 
to the Feynman rules for our diagrammatic expansion. Their presentation will be only 
sketchy and qualitative, since all that will really matter to us in the discussion below is 
counting powers of In fact, we are about to sum an infinite set of leading diagrams 
without ever computing any specific Feynman diagram. 

We start with the propagators. The bare quark propagators (uiU*), {uid*^, (daU*) 
and (dadp) , are simply the four NxN blocks of j^Go- They are all of order j^. (Note that 
the object traced over in (fTBj) and (fT7|) is essentially j^Q^v This trivial observation gives 
at least a practical justification for the factor in the bare quark propagator. A deeper 
reason for this normalization of the quark propagator will be provided below.) These 
propagators are indicated in Fig.l (b)-(e). The diagonal blocks of ^^^o correspond to 
quark propagation without flavor changing. The uu propagator is indicated in Fig. 1(b) 
by a full line, and the dd propagator is indicated in Fig. 1(c) by a dashed line. The 
off-diagonal blocks correspond to quark propagation with flavor changing. The lines are 
all directed since the quarks are "charged" complex flelds. Since all the blocks of j^Qq 
are proportional to the NxN unit matrix, color is preserved under quark propagation, 
whether flavor is changed or not. 




(28) 





(29) 



Non-Hermitean RMT: Planar Diagrams and the "Single-Ring" Theorem 11 

i — ^ j 

(a) 



a a 



(b) (c) (d) (e) 

Fig. 1: The bare gluon propagator (a), and the four blocks of the bare quark propagator (b)-(e). 
All propagators are of order . 



Consider now the "pure glue" action 

Mr1/(0V) = A^m^tr^V + Ar^4tr(0V)^ + NgQtT{(f)^(j)f + ■■■ (30) 

in the exponent of ()20|) . By definition, it is proportional to N, and the "gluon self- 
couplings" m^, g4,gQ,... are all of order A^". Thus, a typical matrix element in 0, drawn 
from the ensemble defined by fj20|) . is of order (and arbitrary phase), which leads, 

by standard arguments, to the conclusion that typical eigenvalues of (p are of order N^, 
as was already mentioned in the context of Ginibre's ensemble. 

The bare gluon propagator can be read-off the quadratic term in (j3n|l in the usual 
manner. We find 

{(pia(p*p)o = J^^V^^P ■ (31) 

It is of order as a direct consequence of the explicit dependence in fpUj) . In order 
that the quark propagator have the same A^-dependence as the gluon's, we scaled the 
quark propagator to be of order as well. This behavior, together with the fact that 
the "gluon self-couplings" m'^,g4,gG, ■ ■ ■ are all of order A^'^, will ensure that the model 
has a smooth behavior as A^ ^ oo, as will be seen in more detail below. 

The bare gluon propagator is depicted in Fig. 1(a). Following 't Hooft [20], it is 
represented by a double line, simply because carries two indices. The fiows along 
these lines preserve both color and fiavor, as is evident from (13 Ij) . The reason why 
the lines in Fig. 1(a) run in opposite directions is easy to understand. Observe that (j) 
transforms under gauge transformations as u ^ S. Thus, the bare propagation of is 
similar to a process in one end of which a w-quark is destroyed and a d-quark is created, 
and in the other end of which the opposite occurs. 

Let us discuss now the vertices. From the "quark-quark-gluon" interaction (|28|) 
we see that a gluon always absorbs a d-quark and emits a w-quark, and (p^ does the 
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opposite. These two processes are described by the two quark-gluon vertices depicted 
in Fig.2(a),(b). Both vertices are taken to be proportional to N, so that they be of 
the same order of magnitude as the gluon self-interaction vertices, such as the quartic 
vertex with coupling Ng4 shown in Fig. 2(c) and the sextic vertex with coupling Ng^ 
in Fig. 2(d). Another way of seeing that the quark-gluon vertices should be of order 
is to note from that with the normalization of the quark propagator we have 



a 



a 



a 



(a) 



(b) 




(c) (d) 

Fig. 2: The two bare quark-gluon vertices (a),(b) and quartic- (c) and sextic- (d) gluon self- 
interaction vertices. All vertices are of order N. 



That all vertices scale in the same way with A^, as do the propagators, is needed to 
obtain a smooth large- A^ behavior. Note from Figs. 2 (a)- (d) that both color and flavor 
flow unobstructed through the various vertices. This is, of course, because all interaction 
terms are gauge invariant. 

To summarize, all bare vertices are of order A^, and all propagators are of order 

Now, use the vertices and propagators defined above as building blocks of graphs - 
Feynman diagrams. The diagrams thus constructed are called double-line diagrams. 

In order to compute a certain observable, all double-line diagrams consistent 
with that observable should be drawn, computed and added with the appropriate 
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combinatoric factors. Note that when a closed index loop occurs in a diagram, its 
index should be summed over. Therefore, we should assign a factor N to each index 
loop. 

We would like now to study the behavior of these double-line diagrams, and show 

that it is determined by their topological properties [201 ■ (An excellent discussion of this 

issue is given in which I shall now follow with the appropriate modifications.) Let 

us start, for simplicity, with connected "vacuum to vacuum" diagrams, i.e., connected 

double-line diagrams without external legs. (Such diagrams arise, for example, in the 

perturbative computation of the normalization factor Z in (j^Uj) . Their sum is essentially 

log Z. The expansion ()27|) of the quark propagator is of course not the sum of vacuum 

diagrams.) Since there are no external legs, each index line in the diagram has to be 

a closed oriented loop. Obviously, there are two types of such loops - loops made of a 

full line and loops made of a dashed line. Let us attach a little surface - a polygon - 

to each of these index loops. Since there are two kinds of loops - there are two kinds 

of such polygons. Thus, we can think of the double-line diagram in question as a tiled, 

triangulated (or more precisely - polygonated) big two-dimensional surface, where any 

pair of edges of two different polygons, which share a gluonic double-line propagator, are 

identified. Since these two edges always run in opposite directions (recall Fig. 1(a)), the 

orientations of all polygons are the same, and the surface thus constructed is orientable jj. 

Moreover, since the two lines which make a gluon propagator are of different types, and 

since index lines flow unobstructed through the vertices, any two neighboring polygons 

which share a gluonic edge must be of different types. In other words, our surface has 

a checkerboard structureff. This surface may have edges, formed by quark lines. Also, 

it will generally be multiply connected - i.e., contain "handles" and internal holes. An 

example of such a (simply connected) surface appears in Fig.3. 

This is all due to the fact that all lines here are directed. For real asymmetric matrices, we would 
obtain unoriented lines and therefore unoriented surfaces e.g., a Klein bottle, 
ft Such checkerboard surfaces appear in certain models of two-dimensional quantum gravity. See e.g. 
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□ 



□ 



MM 



Fig. 3: An example of a simply connected checkerboard surface which contains, in addition to 
quark-gluon vertices, only quartic gluon self-interaction vertices. There is a single quark loop which 
forms the edge. This diagram has 16 quark-gluon vertices plus 16 gluon quartic vertices (V = 32), 25 
index loops (faces, F — 25), and 40 gluon propagators and 16 quark propagators (edges, E = 56). As 
explained below (see (|32|) ). this diagram scales like N. 



Counting the power of N associated with such surface is easy. Assume the surface 
in question has V vertices, E edges and F faces. Each vertex is an interaction vertex 
in the Feynman diagram, and thus contributes a factor A^. Each edge is a propagator 
which carries a factor j^, and each face is an index-loop carrying a factor A^. All-in-all, 
the graph is proportional to ]Sf^~^+^ = N^^ where 

X = F-E + V (32) 

is the Euler character of the graph - a topological invariant of the graph. 

As is well-known, every connected oriented two-dimensional surface is topologically 
equivalent to a sphere with an appropriate number of punctures (holes) and handles 
(or "wormholes") attached to it. For example, a torus is equivalent to a sphere with 
one handle attached, a disk - to a sphere with one puncture, a cylinder - to a doubly 
punctured sphere, etc. Let us denote the number of handles (the so-called genus of 
the graph) by G, and the number of punctures (holes or boundaries) by b. Then a 
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well-known theorem states that 



X = 2-2G-6. 



(33) 



Thus, the leading vacuum diagrams in the limit N oo must have G = b = 0, i.e., the 
topology of a sphere. They are thus proportional to N'^. Obviously, they cannot contain 
quark loops, which are unpaired index loops and thus correspond to holes in the surface 
(6 7^ 0). The last observation is more or less straightforward - one needs just one vertex 
factor to create either a pair of quarks or a triplet of gluons, but there are essentially N 
times more gluon types to sum over than quarks. 

By removing a randomly chosen face off our leading order graph, we can 
continuously stretch it open and project it onto the plane. For this reason, these leading 
vacuum graphs with spherical topology are also known as planar graphs. (We can also 
patch this planar graph back into a sphere, by including the point at infinity.) 

Vacuum diagrams which depend of properties of quarks (i.e., on rj, z, z* in our case) 
must contain at least one quark loop. It follows from (jH^ that the leading such graphs 
in the large-A^ limit must contain precisely one quark loop (i.e. - a single hole) and no 
handles. Such graphs have the topology of a disk, and could thus be drawn on the plane. 
They are planar graphs, like the pure glue graphs, with the sole difference between the 
two cases being that the outer boundary of that graph is that of the hole - i.e., the quark 
loop. There can be no quark loops in the interior of the graph in the large- limit. 

The reader may identify the Feynman rules described above as those which 
correspond to the zero- dimensional "quantum field theory" with lagrangian 



Diagrams which contribute in higher orders of the expansion of this theory do contain 
quark loops in their interiors, i.e., correspond to surfaces with internal holes. Diagrams 
with such internal quark loops obviously do not appear in the expansion of j^i^Q^j,) since 
the random matrix probability distribution ()20|) really does not contain any dynamical 
quarks. One method of suppressing these internal quark loops to any order in the 
expansion of (jMj) is to apply the so-called "replica trick" f Since here we are 

interested only in the N oo limit, namely, the leading order in the expansion, there 
is no need to invoke any replica considerations. 

f That is to say, to replicate the quarks in H34|l into n copies, with the same inverse bare propagator and 
couphng to the gluons. One then computes various observables in the replicated theory as a function 
of n, considered as a complex variable, and then take the limit n ^ 0. Each internal quark loop costs 
a factor n and is therefore suppressed in that limit. 




(34) 
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3.1.2. Computation of j^Q^v to Leading Order in Obviously, by cutting open the 
quark loop at the boundary of our vacuum diagram at one point, and attaching color 
indices /i and u at the two ends, we obtain a diagram which contributes to j^G^u in the 
expansion (j2Z|)- Thus, to leading order in j^, all diagrams which contribute to j^G^u 
(and are therefore of order j^) are planar diagrams, with the quark line running at 
the boundary, and all gluon double-lines are attached to it on the same side, without 
crossing each other. A couple of such diagrams appear in Fig. 4. 





Fig. 4: Two simple planar diagrams which contribute to j^Gfn^- The upper one contributes to 
the uu block and the lower one to the dd block. Both diagrams are of order Color indices and 
arrows are suppressed. These diagrams are one-quark-reducible: their topology is such that they can 
be split into two disjoint pieces by cutting a single quark line. 



I shall now explain how to sum all these planar diagrams systematically and thus 
obtain j^Gf^u in the limit iV — > oo. 

The quark propagator p4|) is given by the Dyson- Schwinger equation in terms of 
the one-quark irreducible self-energy NT,^^: 

We emphasize that this equation always holds, and practically amounts to the definition 
of the one-quark irreducible self energy. Evidently, all planar self-energy diagrams are 
of order N\. Thus, to leading order, S^,^ is of order A^°. A couple of planar one-quark 
irreducible diagrams which contribute to NH^j^ appear in Fig.5. 

I Their large- iV behavior is like that of the planar diagrams which contribute to j^G^u with their 
two external legs amputated. Amputating two external legs amounts to multiplication of the latter 
by A''^. Another way of looking at this is to say that the self-energy should scale with A^ as does the 
quark-quark-gluon vertex NV which is combined with the inverse propagator in the averaged quantity 
in imi). 
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Fig. 5: Two simple planar one-quark irreducible diagrams which contribute to the uu block of 

NTifii^: a diagram with no gluon self-interactions (a), and a diagram with a quartic gluon self-interaction 
(b). Evidently, both diagrams are proportional to 6ij. Also, both diagrams are of order N. 



For pedagogical clarity, we also display in Fig. 6 a couple of non-planar self-energy 
diagrams which contribute in the next-to-leading order. 




(a) (b) 

Fig. 6: Two simple non-planar self-energy diagrams with a single crossing of gluon lines. Color 
indices and arrows are suppressed. The effect of this crossing is dramatic: in diagram (a) there are no 
index loops at all (as opposed to its (reducible) counterpart in the upper part of Fig. 4 with two such 
loops), and in (b) there is only one index loop (as opposed to its planar counterpart in Fig. 5(b) with 
three such loops). These diagrams are both of order i.e., down by a factor iV~^ with respect to the 
leading order. 



The planar self-energy diagrams shown on Fig. 5 are clearly proportional to the unit 
matrix. By carefully studying a few additional planar self-energy diagrams the reader 
can convince himself or herself that this is a general property: in these planar self- 
energy diagrams, there is always a unique peripheral index line, which emanates from 
the quark-quark-gluon vertex at one end of the diagram, and participates in several 
gluon propagators and flows unobstructed through a number of gluon self-interaction 
vertices, without crossing any other index line, all the way to the quark-quark-gluon 
vertex at the other end of the diagram. (A slightly different way of looking at this is to 
connect the two ends of the diagram by a quark propagator and close it into a planar 
vacuum diagram. The point of gluing lies on a quark line, which is of course the outer 
edge of one of the polygons making the surface. The other edges of that polygon lie 
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within the surface and combine to a directed index hne, which becomes the peripheral 
index hne mentioned above upon cutting the diagram open.) Thus, evidently, the off- 
diagonal blocks of NTi^y are null and its uu and dd blocks must be proportional to 
the unit matrix. Furthermore, the u and d quarks are completely equivalent. (We can 
switch the roles of (j) and 0^ in ()20|) and nothing would change.) Thus, the uu and dd 
blocks of Tj^i, must be equal. In other words, for probability distributions of the form 
(EOl) 

S^j, = S(?7;r)(5^^ (36) 

(where we recall that r = \z\.) Note that S depends on r and not separately on z and 
z* because it may be expanded as an infinite sum of traces of Qq, each term of which is 
a function only of rj and r. Thus, from ()35p . 

( Tj — ri z \ 



r"^ — {rj — S)" 



(37) 



\^ z* ^-v J 



leading to 



^ = ^E^M.-^2_^( _\)2- (38) 

The important point |3I211 is that in the large- limit, the one-quark irreducible self- 
energy NT, can be written in terms of the cumulants of P{4') = (1/Z)e~^*''^*^'^^'^^ 
(Eq. ()20|l ). namely, the connected correlators involving k 0's and 0^'s. Thus, to leading 
order in j^, this cumulant is the sum of all connected planar pure glue diagrams with 
k external 0's and 0^'s (with a prescribed index structure). By definition, is a 
quantity of order A^. (It could be thought of essentially as the "dressed" gluon self- 
interaction vertices Ng2k-) This makes sense, since any of the diagrams contributing to 
can be thought of as a planar, connected pure glue vacuum diagram (and thus of 
order A^^), cut open at k places along its peripheral gluon propagators§ in such a way 
that it remains connected (i.e., there must be at least one gluon self-interaction vertex 
between consecutive cutting points), and then amputating the gluon legs thus formed. 
A moment's reflections reveals that this operation removes k + 1 index loops (i.e., k 
internal index loops as well as the outer perimeter loop) and also k gluon propagators. 
Thus, it costs a factor of j^, reducing the A^-power counting of the diagram from 2 to 1. 

§ These cuts must be made along peripheral gluon propagators since the gluon diagram in question 
must be connected, without crossings, to other parts of a bigger planar diagram. 
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Diagrammatically speaking, the may be thought of as a proper vertex of k 0's 
and 0^'s, i.e., upon embedding it in a planar diagram it appears as a "blob" out of which 
emanate k 0's and 0^'s in an alternate manner and without crossings, and which cannot 
be separated into two smaller blobs, with ki 0's and 0^'s in one blob and /c2 0's and ^''''s 
in the other (with ki + k2 = k of course.) 

The way NT, is determined by the cumulants can be deduced by working out 
a few simple cases. For example, the contribution of r2 is shown in Fig. 7. 



Fig. 7 : Contribution of the quadratic cumulant r2 to self-energy. 

Counting powers of in Fig. 7 is as follows: r2 is of order A^. The factors 
associated with the two gluon propagators poking out of the dark gluon blob representing 
r2 are canceled against the factors of in the quark-gluon vertices. The horizontal line 
at the bottom is the dressed planar quark propagator, i.e., the sum of all planar diagrams 
contributing to j^Q^u, and is of order j^. It is connected to the internal index line which 
emanates from the quark-quark-gluon vertex on the left and flows through the gluon 
dark blob. On the way, that line participates in several gluon propagators and flows 
unobstructed through a number of gluon self-interaction vertices in any of the planar 
pure glue diagrams which make the blob, without crossing any other index line, all the 
way to the quark-quark-gluon vertex on the right. Thus, it preserves its color index, 
which must be summed over. Thus, the dressed quark propagator is traced over and 
produces a factor Q{ri;r) of order N^, as defined in (|38|) . The overall contribution of 
Fig. 7 to self-energy is therefore r2^ which is of order A^. Evidently, the diagram shown 




J 
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in Fig.5(a) is the leading contribution to this quantity. It will be instructive to analyze 
the contribution of the next cumulant, as well. Its contribution to self-energy is 
shown in Fig. 8. 




Fig. 8: Contribution of the quartic cumulant r4 to self-energy. 



Counting powers of in Fig. 8 is as follows: r4 is of order A^. The factors 
associated with the four gluon propagators poking out of the dark gluon blob are 
canceled against the factors of in the quark-gluon vertices. Each one of the three 
dressed quark propagators at the bottom is connected to an internal index loop which 
flows through the gluon dark blob. Thus it is traced over and produces a factor Q{ri] r). 
The overall contribution of Fig. 8 to self-energy is therefore T4^Q^ which is of order A^. 
Evidently, the diagram shown in Fig.5(b) is the leading contribution to this quantity. 

The general pattern should be obvious||. Thus, we obtain [2J 

oo 

NE{r);r) = Y,^,,[g{r^;r)f-\ (39) 

k=l 

where each of the summed terms is manifestly of order A^. 

The requirement, mentioned above, that the blob corresponding to r2A: be 

inseparable into smaller blobs (a property referred to as "gluon connectedness" in [21]) 

can be explained with the help of the simple example shown in Fig. 9. 

II It should also be clear that we need not dress the quark-quark-giuon vertex, as such corrections are 
already accounted for by the terms summed in (|39() . 




fa) fb) 

Fig. 9: The contribution of the quadratic cumulant (a) and that of a separable gluon blob (b). 
The diagram (b) is evidently already accounted for in (a). 



Evidently, the diagram in Fig. 9(b) is already accounted for by the contribution 
of the second cumulant. It is therefore not included in the contribution of the fourth 
cumulant shown in Fig. 8. 

If we knew r) in terms of Q{r]] r), then we can solve for r) and hence the 
density of eigenvalues. We are thus faced with the problem of determining the r2fc's 
which appears to be a rather difficult task^. 

It is important to note that the in flH^ depend only on the probability 
distribution P{(f)) and not on the particular quantity we average over. Consider the 
problem of determining the eigenvalue density of the hermitean matrix We would 
study the Green's function 

oo 

F{w) = — tr,^, — = / ^ — , 40 

^ ^ W w-(l)^' J w-a ^ ^ 



where is the averaged eigenvalue density of We see from p7p that this Green's 
function is related to Q{ri;r = 0) simply through 

g(n;r = 0) =riF{ri'^). (41) 

Also, for z = we have from (j38p 

nv;0)=v-G~\v;0). (42) 

The crucial observation is that F{w) is already known in the literature on chiral and 
rectangular block random hermitean matrices for the Gaussian distribution, as well as 
for non-Gaussian probability distributions of the form (j2(jp with an arbitrary polynomial 
% Except, of course, for the Gaussian ensemble, for which only the variance r2 ^ 0, and thus simply 

s = r2g. 
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potential V^(0V))[I3IIHllini- (See also the paper by Feinberg and Zee cited in [ini-) In 
fact, the authors of [IHl HH] simply calculated the diagonal elements of the propagator 
OfMuiVi^) using Dyson gas techniques (so the coefficients r2fc are only implicit in these 
papers.) 

The contents of (jH^ . (PT|) and (jl^ may thus be summarized as follows: there is a 
unique function a(^) which behaves like 1/^ as ^ tends to infinity, is regular at ^ = 0, 
and satisfies the equation ^ — = YlT=i ^2kO''^^~^- This function is 

a(^)=^Fie)^g{^;0). (43) 

We, on the other hand, are interested in the opposite case, where rj = 0,r = \z\ 7^ 0. 
Consider the matrix ^^jy(0 — iO; z, z*). It follows from and that its x blocks 
are all proportional to the unit matrix l^r (we refer to this as "color index democracy"), 
and we see that 

g G*\ 



g^,„iO - iO;z,z*) 



(44) 



G g J 

where ^(0 — iO; r) is given in (jl8p and the complex Green's function G{z, z*), the lower 
left block of g,_iu{0 — ^0; z, z*), was defined + in (jH)). Inverting the propagator in (ji^ and 
using (fT^ (with rj = 0) and ((SHI), we obtain the two equations 

G* 



(45) 

|Ur|- — y- 

and 



\G\^-g^ 
g 



Note from ()45|1 that if ^ = 0, then G = 1/z, which corresponds to the region outside 
the domain of the eigenvalue distribution. In addition, using twice, we obtain 

zG = = {\Gf - g') r\ (47) 

Thus, due to the fact that g is pure imaginary (see (fTHjl ) we conclude that 7 = zG is 
always real and non-negative. Alternatively, we can write 

(48) 

and therefore from the fact that < we conclude that < 7 < 1, in accordance with 

the monotonicity of 7(r) and the sum-rules 7(0) = and 'y{oo) = 1 in ()24|) . 

+ From this point on, unless otherwise stated, our notation wiU be such that Q = 5(0— iO; r), E = S(0; r) 
and G = G{z, z*). Note in particular from (|18|l that Q is then pure imaginary, with "^Q > 0. 
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(49) 



k=l 



where the last equahty follows from (jH^ (with rj = 0). Using and PU)) we thus 
have simply 



, zG 



(50) 



It then follows from PHj) and the statement preceding that equation, and from (I^Uj) that 

me) = G = ^- (51) 



We have thus eliminated the ! 

From ()48|) and from ()50|) we determine = 7/^ to be 

e 



7 2 

r 



7 



(52) 



Note that < 0, consistent with the bounds < 7 < 1. Comparing and we 
obtain the remarkable relation 

g{^;0) = g{0-r), (53) 

for ^ and r that are related by ()52j) (with ^ imaginary.) 

We now have the desired equation for 7 = zG: substituting into we obtain 



7 



^2 ^ 



7 r 
7 - ■ 



7 + 1 



0. 



(54) 



This is the master equation announced in the Abstract. Thus, given F we can solve for 
7(r) using this equation. 

From ()40|) we observe that for w large we may expand 

F(w) = - + Xf(w), (55) 



w 



where f{w) = fo + ifi/w) + ■ ■ ■. Substituting w = = jr"^/ (7 — 1) into and letting 
7 tend to 1, becomes 

7H 



(7-1)^ 



7r^ 



(56) 



and thus 7 = 1, which corresponds to the region outside the eigenvalue distribution, is 
always a double root of (jSH). 

We will now explain how the Green's function F{w) is obtained in the literature |18| 
HH]. (See also the paper by Feinberg and Zee cited in [IE]-) Let V he a. polynomial of 
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degree p. From the definition of F{w) in (j4Uj) and its analyticity property, we expect F, 
following the arguments of Brezin et al, to have the form* 

F{w) = ^V'{w) - P{w) ^{w-a){w-b) , (57) 

where 

p-2 

P{w)=Y,<^kwK (58) 

k=-l 

The constants < a < 6 and Ck are then determined completely by the requirement 
that F{w) ^ ^ as 1(7 tends to infinity, and by the condition that F{w) has at most an 
integrable singularity as w — > 0. (Thus, if a > 0, inevitably c_i = 0. However, if a = 0, 
then c_i will be determined by the first condition.) 

4. An Example: the Quartic Ensemble 

Having determined F[w) in this way, we substitute it into ()54|) and find G{z,z*). 
Equation ()54|) is an algebraic equation for 7(r) and thus may have several r dependent 
solutions. In constructing the actual 7(r) one may have to match these solutions 
smoothly into a single function which increases monotonically from 7(0) = to 
7(00) = 1 (recall ()24|) ). (An explicit non-trivial example of such a sewing procedure 
is the construction of 7(r) in the disk phase of the quartic ensemble[2].) We can thus 
calculate the density of eigenvalues explicitly for an arbitrary V . 

The Gaussian case with V = (f)^(f) provides the simplest example. Here we have 
2^F(^2) ^ ^ _ ^^2 _ 4^ whence the roots of §^ are 7 = 0, 1 and r^. We note that 
7 = is unphysical, 7 = {i.e., G = z*) corresponds to Ginibre's disk ^14], and 7 = 1 
is the solution outside the disk. 

As an explicit non-trivial example, the quartic ensemble with V{(j)'^(f)) = 2m^0V + 
(7(0^)^ studied in detail in [2]. As should perhaps be expected in advance, the 
following behavior in the parameter space m'^,g > was found: for positive, the 
eigenvalue distribution was disk-like (and non-uniform), generalizing Ginibre's work, but 
as = — /i^ was made more and more negative, a phase transition at the critical value 
fil = \/2g occurred, after which the disk fragmented into an annulus. The density of 
eigenvalues was calculated in in detail. In particular, the radii of the disk and annulus 
were calculated explicitly in terms of the couplings in V . An interesting behavior of the 
density of eigenvalues was observed through the disk-annulus phase transition. As the 
* Here we assume for simplicity that the eigenvalues of (/)^(/) condense into a single segment [a, 5]. 
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transition was approached from the disk side (e.g., by lowering towards — /i^ holding 
g fixed), a central region of finite radius in the disk got progressively depleted. At the 
transition depletion was complete and an annulus of finite critical radius was formed, 
with vanishing eigenvalue density at the inner radius of the annulus. As the couplings 
were tuned deeper into the annular phase, a finite jump in the density of eigenvalues 
at the inner radius appeared. The interested reader could find all relevant formulas in 
Section 5 of 2j and Section 4 of jS]. In what follows we shall compare these analytical 
large- results, following [3], against numerical simulations. 

In Fig. 10 we display the numerical results for pdisk{f) for 128x128 dimensional 
matrices, and compare them to the analytical large- A^ result (5.8) of 2J (or (4.4) of [H]). 
(As a trivial check of our numerical code, we also included in this figure the results for 
the gaussian (Ginibre) ensemble.) Evidently, the agreement between the numerical and 
the analytical results is good. Note the finite- A^ effects near the edge of the disk. 



N=128 



^ I I I I I I I I r 




r 

Fig. 10: Comparison between Monte-Carlo measurements of the density of eigenvalues p{r) of 
matrices (j) of size 128x128, taken from the quartic ensemble V{(f>''(j)) = 2m'^(j)''(j) + g{(l)''(j))'^ with = 0.5 
(disk phase) and for g — 0,0.5,1,2,4 {g increases from bottom to top), compared to the analytical 
results of ^ (solid lines). At 5 = we obtain Ginibre's Gaussian ensemble with V = 4>^4>, with its unit 
disk of eigenvalues. (This figure is taken from (3|.) 
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In Figsll.(a)-(c) we display our numerical results for Pannuiusij) for matrices of 
various sizes, and compare them to the analytical large- iV result (5.18) of j2] (or (4.13) 
of [3]). In these figures we hold —m? = fi^ = 0.5 fixed, and increase g from 0.025 to 0.1. 
(Here we have /x^ = 0.5 = fil/2y/2g. Thus increasing g as indicated in the text brings 
us closer to the disk-annulus phase transition.) 

N = 64.1 28,256 



lo in in 




4 4.2 4.4 4.6 4.8 2.6 2.8 3 3.2 3.4 3.6 1.6 1.8 2 2.2 2.4 

r r r 

Fig. 11: Results of Monte-Carlo measurements of the density of eigenvalues p(r) of matrices (j) of 
sizes corresponding to iV = 64, 128 and 256, taken from the quartic ensemble with = — /j,^ = —0.5 
(annular phase) for various values of the quartic couphng: g = 0.025 in [a], g = 0.05 in [b] and g — 0.1 
in [c]. These are compared to the analytical results of [21 (solid lines). As N increases, the numerical 
results converge monotonically to the analytical results. (This figure is taken from [H].) 

In jH] we have also studied numerically the disk-annulus phase transition. We have 
measured the density of eigenvalues p(r) of matrices of size 128x128, taken from the 
the quartic ensemble with /x^ = 0.5 and for g = 0.025,0.05,0.1,0.125,0.15 and 0.175. 
The results are displayed on Figure 12. 

For these values of g, we start in the annular phase at the lowest value of g. For 
our set of parameters we have fi'^ = 0.5 = ^^I2^j2g. Thus, increasing g (while keeping 
[]? fixed at 0.5) brings us closer to the disk-annulus phase transition, which occurs (at 
large A^) at gc = 0.125. Increasing g beyond that, puts us into the disk phase. 

The first three profiles on the right in Figure 12 belong to the annular phase. As g 
increases towards the transition point at gc = 0.125, these three graphs exhibit decrease 
of the inner radius of the annulus, in accordance with the analytical results. 

The critical density profile, corresponding to gc = 0.125, is the fourth profile (from 
the right). For our choice of parameters, the theoretical boundary radii of the critical 
annulus, i.e., at g = 0.125, are i?f^** = l//ic = and = V^/a^c = 2. These 
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boundary values fit nicely with the features of the critical profile in Figure 12. 

Finally, the last two profiles in Figure 12 have pronounced tails extending to r = 
and thus belong to the disk phase. 




Fig. 12: Monte-Carlo measurements of the density of eigenvalues p{r) of matrices (j> of size 
128x128, taken from the the quartic ensemble with fi"^ = 0.5 and for g = 0.025, 0.05, 0.1, 0.125, 0.15 and 
0.175 {g increases from right to left). The first three profiles on the right (corresponding to the three 
lowest values of g) evidently belong to the annular phase. The fourth density profile from the right is 
the critical one (corresponding to gc = 0.125). Finally, The last two profiles (which correspond to the 
two higher values of g) belong to the disk configuration. (This figure is taken from [5].) 



5. Boundaries and Shape Universality of the Eigenvalue Distribution: the 
Single Ring Theorem 

A remarkable property of ()54|) is that it has only two r-independent solutions: 7 = 
and 7 = 1 12]. Since the actual 7(r) increases monotonically from 7(0) = to 7(00) = 1, 
without any plateaus in between, we immediately conclude from this observation that 
there can be no more than a single void in the eigenvalue distribution. Thus, in the 
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class of models governed by P{4>) = le"^*''^^''^^'^) (Eq. ((201)), the shape of the eigenvalue 
distribution is either a disk or an annulus, whatever polynomial the potential V{(f)^(f)) 
is. This result is the "Single Ring Theorem" of |2J. 

The "Single Ring Theorem" may appear counter-intuitive at first sight. Indeed, 
consider a potential V{(f)'^(f)) with several wells or minima. For deep enough wells, we 
expect the eigenvalues of 0V "fall into the wells". Thus, one might suppose that 
the eigenvalue distribution of (p to be bounded by a set of concentric circles of radii 
< Ti < r2 < ■ ■ ■ < Tn^^^, separating annular regions on which p(r) > from voids 
(annuli in which p(r) = 0.) A priori, it is natural to assume that the maximal number 
of such circular boundaries should grow with the degree of V, because V may then have 
many deep minima. Remarkably, however, according to the "Single Ring Theorem" the 
number of these boundaries is two at the most. 

To reconcile this conclusion with the a priori expectation just mentioned, note 
that while the eigenvalues of the hermitean matrix may split into several disjoint 
segments along the positive real axis, this does not necessarily constrain the eigenvalues 
of itself to condense into annuli. Indeed, the hermitean matrix (f)^(j) can always be 
diagonalized (p'^cf) = WA'^U by a unitary matrix U, with = diag(A^, A2, ■ ■ ■ , A^), 
where the Aj are all real. This implies that = V'^AU, with V a unitary matrix as well. 
Thus, the complex eigenvalues of (p are given by the roots of det(2; — AW) = 0, with 
W = UV'^ . Evidently, as W ranges over U{N) (which is what we expect to happen in 
the generic case), the eigenvalues of AW could be smeared (in the sense that they would 
not span narrow annuli around the circles \z\ = \Xi\.) 

The last argument in favor of the "Single Ring Theorem" clearly breaks down 
when W fails to range over U{N), which occurs when the unitary matrices U and V 
are correlated. For example, may be such that W = UV"^ is block diagonal, with the 
upper diagonal block being a K x K unitary diagonal matrix diag(e*'^^, ■ ■ ■ , e*'^^) (and 
with K a finite fraction of A^). In the extreme case A' = A^, in which W is completely 
diagonal, W = e^'^ = diag(e*'^^, ■ ■ ■ , e*'^^), we see that = We'^'^AU is a normal matrix 
(that is, [0,0^] = 0) with eigenvalues diag(e*'^iAi, ■ ■ ■ ,e*'^^A7v)- Thus, normal matrices, 
or partially normal matrices (i.e., the case K < N), evade the "Single Ring" theorem: if 
the first K eigenvalues A^, A|, ■ ■ ■ , A|^ of 0^ split into several disjoint segments along the 
positive real axis, the corresponding eigenvalues of will split into concentric annuli in 
the complex plane obtained by revolving those A-segments. Normal, or partially normal 
matrices are, of course, extremely rare in the ensembles of non-hermitean matrices 
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studied in this paper, and do not affect the "Single Ring" behavior of the bulk of 
matrices in the ensemble. 

Let us sketch now the proof of the Single Ring Theorem, following [2]. Thus, let us 
assume for the moment that the domain of eigenvalues has n boundaries. It is easy to 
see from ()22|) . the defining equation of 7(r), that in the annular void Vk < r < r^+i in 
the eigenvalue distribution, 7 is a constant which is equal to the fraction of eigenvalues 
contained inside the disk r < r^. Thus, the equation to determine 7 (Eq. (j^^ which 
we repeat here for convenience) 



7 



r^F(-^)-7 + l 



(59) 



must have a series of monotonically increasing constant solutions 71 < 72 < • ■ • < 1, 
which correspond to the various voids. 

In particular, from (j24|) we have 7 = 1 for r > r^n.^^, namely, G = 1/z. Thus, 7 = 1 
must be a solution of ()59j) . the maximal allowed constant solution, as we already saw 
following Also, if vi > 0, namely, if there is a hole at the center of the eigenvalue 
distribution, then for r < ri we must have 7 = (independently of r), which is obviously 
a solution of (jSHI)- (From the paragraph right below (j37|) it is clear that wF{w) — > 
at w = 0. Thus 7 = is indeed a solution. On the other hand, if ri = 0, so that the 
eigenvalue distribution includes the origin, an r independent solution 7 = is of course 
only a spurious solution which should be discarded.) 

Assume now that 7 = 70 is an r independent solution of (j59j) . Taking the derivative 
of (j^Uj) with respect to at 7 = 70 we obtain 

,dF' 



7o 



= 0, (60) 

70-1 

which is the condition for the existence of an r independent solution 70. Thus, there 
are two possibilities: either 70 = 0, which we already encountered, or F{^'^) +^^^ = 0. 
This equation immediately yields 

Fie) = ^,, (61) 

where the integration constant is fixed by the asymptotic behavior of Fl^"^) as ^ 00. 

But for a generic V{(j)^(j)), F{w) is given by (Eq. (HO}) 

00 

^/ N , 1 1 , f p{(j)da 
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with (jUT|) being the asymptotic behavior of F{w) as w = ^ oo. We thus conclude 
that ^ oo in ()6H). namely, that 70 = 1. Thus, to summarize, the only possible r 
independent solutions of are 7 = and 7 = 1, which we already discussed. Since 
no other r independent solutions arise, there can be no more than a single void in the 
eigenvalue distribution, whatever polynomial the potential V{(j)^(j)) is. The shape of the 
eigenvalue distribution is thus either a disk or an annulus. 

Interestingly enough, we can arrive at the same conclusion by invoking other general 
aspects of the method of hermitization, and thus providing a nice self consistency check 
of our formalism. Recall at this point that the boundaries of the eigenvalue distribution 
are given by (fT^. namely, the zeros of Q. But Q is given by (Eq. (HH)) 
r2 7(7-1) 

y = — ^ — , 

from which we see that at the boundaries 7 can take on only two values : zero and one. 
Since 7 is a constant in the void, by continuity, these are the only possible values of 7 
inside any of the voids. Therefore, there may be two circular boundaries at most. In 
other words, as far as the eigenvalue density of (p is concerned, an ensemble of the form 
(pn)) may have two phases at most, as we concluded earher. 

In addition to 7 = 1 (and possibly 7 = 0), must have other roots which 
do depend on r. Among all these other roots, we expect to find a unique root 7(r), 
which is a positive monotonically increasing function of r, that matches continuously 
at the boundaries {k < 2) to the constant roots of (|59|). The actual zG{z,z*) of 
the ensemble (l2Up is therefore a continuous monotonically increasing function of r, as 
required by (j22j) . If the eigenvalue distribution is annular with boundaries < ri < r2 
(the case vi = corresponding to the full disk), it vanishes for < r < ri, rises 
monotonically from zero to one on ri < r < r2, and equals to one for r > r2. 

The sextic potential l^(0V) = "^^^V + |(0V)^ + f is the potential of lowest 

degree in (|2Up for which the eigenvalues of 0V split into more then a single segment. 
Therefore, it is the simplest case in which the Single Ring theorem has non-trivial 
consequences. In fact, it is easy to see that there can be at most two eigenvalue segments 
in the spectrum of iii which case, one segment touches the origin. In |3j we used 
Monte-Carlo simulations to generate an ensemble of random matrices (of size 32 x 32) 
corresponding to ((201) with the sextic potential. The coefficients of the latter were tuned 
to the phase where the spectrum of splits into two segments. We verified numerically 
that the eigenvalue distribution of in this phase was always a disk, and not a disk 
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surrounded by a concentric annulus. 
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6. Universal Features of the Disk and Annular Phases and the Transition 
Between Them 

According to the "Single Ring" theorem, the eigenvalue distribution of (p is either a disk 
or an annulus. The behavior of F{w) as w ~ turns out to be an indicator as to which 
phase of the two the system is in [3], as 1 now show: 

A. Disk Phase: In the disk phase we expect that p(0) > 0, as in Ginibre's case. 
Thus, from (PHj) p(r) = {l/r){d'~f/dr) = 2{d'~f/dr'^) and from the first sum rule 7(0) = 
in we conclude that 

7(r) ~ ^piOy (62) 

near r = 0. Therefore, for r small, ()54|) yields 

-^- + ---J~-^, (63) 

namely, F{w) ~ f*^^ w ~ 0, as we could have anticipated from Ginibre's case. 

This means that in the disk phase we must set a = in ()57|). Consequently, in the disk 
phase c_i does not vanish. We can do even better: paying attention to the coefficients 
in (j37j) and ()58|) (with a = 0) we immediately obtain from (jUHj) that 



B. Annular Phase: In the annular phase 7(r) must clearly vanish identically in the 
inner void of the annulus. Thus, (j5l|l implies that F{w) cannot have a pole aX w = 0, 
and therefore from we must have c^i^/ab = 0. Thus, the annulus must arise for 
c_i = (the other possible solution a = 0, c_i 7^ leads to a disk configuration with 
7 = only at r = 0, as we just discussed.) 

Thus, to summarize, in the disk phase F has the form 



FdiskM = ^V'iw)- \^^P^w-' + Co + CiW + --- + Cp_2 w^'^j ^/wiw - b) , (65) 
while in the annular phase it has the form 

Fannulus{w) = (w) - (cq + Ci W H h Cp_2 W^'^) a/(w - a){w - h) . (66) 

Having determined F{w) in this way, i.e., having determined the various unknown 
parameters in or in (jMj) . we substitute it into (jMj) and find G{z,z*). We can 
thus calculate the density of eigenvalues p{r) explicitly for an arbitrary V . 
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C. The Disk-Annulus Phase Transition: We now turn to the disk-annulus phase 
transition. An important feature of this transition is that F{w) is continuous through 
it. To see this we argue as follows: By tuning the couplings in V , we can induce 
a phase transition from the disk phase into the annular phase, or vice versa. Note, 
of course, that we can parametrize any point in the disk phase either by the set of 
couplings in V or by the set of parameters {c_i, Cq, ■ ■ ■ Cp_2; &} in ()65|). The "coordinate 
transformation" between these two sets of parameters is encoded in the asymptotic 
behavior of F{w). Similarly, we can parametrize any point in the annular phase either 
by the set of couplings in V or by the set of parameters {cq, ci ■ ■ ■ Cp-2] a, b} in (jMj) . Due 
to the one-to-one relation (in a given phase, once we have established it is the stable 
one) between the couplings in V and the parameters in F{w) — ^V'{w) (namely, the 
c„'s and the locations of the branch points of F{w)), we can describe the disk-annulus 
transition in terms of the latter parameters (instead of the couplings in V). Clearly, the 
transition point is reached from the disk phase when p(0) = 0, that is, when c_i in 
vanishes: 

c^rf = . (67) 

Similarly, the transition point is reached from the annular phase when the lower branch 
point a in ()66|) vanishes. Thus, e.g., in a transition from the disk phase into the annular 
phase, Fdisk{w) in would cross-over continuously into Fannuiusiw) in (jUBj) through a 
critical form 

The continuity of F{w) through the transition was demonstrated explicitly in j2] 
for the quartic ensemble V{(f)^(f)) = 2m^0^(/) + g{(j)'^(j))'^. 

This discussion obviously generalizes to cases when F{w) has multiple cuts, which 
correspond to condensation of the eigenvalues of into many segments. If w = is a. 
branch point of F{w), that is, if the lowest cut extends to the origin, we are in the disk 
phase, 

Fdisk{w) = ^V'{w)-{c^iw'^ + Co + ci w H h Cp_2 w^"^) a/w(w - 61) ■ ■ ■ (u; - 6„) , 

with < bi < ■ ■ ■ < bn- The relation (jMj) then generalizes to 

/ P(0) 

Since c_i must be real we conclude that such a configuration exists only for n odd. 
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If the lowest branch point in F{w) is positive, we are in the annular phase with 

Fannulus{w) = ^V\w)-{co + Ci W ^ h Cp_2 W^'^) a/(w - o) (w - 61 ) ■ ■ ■ (w - 6„) 

The phase transition would occur when the couplings in V{(j)''(f)) are tuned such 
that Fdiskiw) and Fannuiusiw) match continuously, as was described in the previous 
paragraph. 

7. Boundaries and Boundary Values 

Remarkably, with a minimal amount of effort, and based on the mere definition of F{w) 
(Eq. (j4Uj) ) ■ we were able to derive in jH] simple expressions for the location of the 
boundaries of the eigenvalue distribution and also for the boundary values of p(r) in 
terms of the moments of p(cr), which, I remind the reader, is the density of eigenvalues 
for a hermitean matrix problem. 

To this end it is useful to rewrite our master formula (j54p for 7(r) as 

wF{w) = 7 (72) 

with 

w = . 73 

7 - 1 

We start with the outer edge r = Rout (either in the disk phase or in the annular phase.) 
Near the outer edge 7 — > 1—, and thus w —00. We therefore expand F{w) in (f72|) in 
powers of 1/w. As can be seen from fj4Up . the coefficient of w~^~^ in this expansion is 
the k-th moment 

00 

<a'' >= I p{a)a''da (74) 





of p((j) (which is of course normalized to 1.) For the class of models we are interested 
in here, all the moments < >, k > are finite. From the expanded form of (j72|) one 
can obtain, by a straightforward calculation that 

RLt =<^> (75) 

and 



p{Rout) = . (76) 

< > _ < o- >2 

Thus, is simply the first moment of p(cr), and the density p{Rout) is inversely 

proportional to the variance of a. 
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For the p{a) under consideration here, < o"^ >, and consequently p{Rout), are always 
finite. Outside the boundary p(r) vanishes identically, of course, and thus, p(r) always 
"falls off a cliff" at the boundary, for all probability distributions of the form (pUj) with V 
polynomial. It would be thus interesting to study circularly invariant matrix ensembles 
P{(j)'^(f)) such that the eigenvalue distribution p(cr) of 0V has a finite < a > but an 
infinite < cr^ >. Then p{Rout) would vanish. This would naturally raise the question 
whether in such situations, p(r) behaves universally near the edge (that is, if near the 
edge it vanishes like {Rout ~ f^Y with e being some universal exponent). 

We now turn to the annular phase, and focus on the inner edge r = i?j„ of the 
annulus. According to the discussion in the previous section (see Eq. ()66p and the 
discussion above it), a > in (jSZj), and thus F{w) is analytic in the domain \w\ < a. 
We therefore expand F{w) in ()72|) in powers of w. As can be seen from ()40|) . the 
coefficient of in this expansion is the negative moment 

oo 

i—k) = / Pannulus{(y)—rda (77) 



a J a 



of p annulus {o') . From (jFTf) (or (jU^ ) it is possible to see that all these negative moments 
of Pannuius{cr) cxist. From this expanded form of (f?^ one can obtain, again, by a 
straightforward calculation [3^, that 

(-) (78) 



2 

"in 

and 



= ^ -1^2 ■ (79) 

< a > — < a ^ 

Thus, R^^ is simply the moment of Pannuiusic) , and the density p{Rout) is inversely 
proportional to the variance of cr^^ . Since <l/cr^>, k = 1,2 are finite, Pannuius{cr) 
jumps from zero (in the inner void of the annulus) to a finite value at the inner edge 
Rin- It can be shown however, that when a — 0, that is, in the annular to disk 
transition, < 1/a > remains finite, but < l/a^ > diverges like Xj^fa. Thus, from (fTSjl 
we see that Rinnerio- = 0), the critical inner radius, is finite. The annulus starts up 
with a finite inner radius. Also, in this limit, we see from (fTUj) that p{Rin) vanishes like 
^/a. As we approach the annulus-disk transition, the discontinuity in p(r) at the (finite) 
inner edge disappears. 

We saw at the end of the previous section (see Eqs. (jHSj) -(jHH)) that F{w) 
is continuous through the disk-annulus phase transition. Thus, our master formula 
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wF{w) = 7 to determine 7(r) (Eq. (f72|)) is also continuous through the transition. 
Consequently, p(r) = {l/r){d'y/dr) must remain continuous through the disk-annulus 
transition, and has (at the transition) the universal behavior described in the previous 
paragraph. 
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